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The escape probability out of the superconducting state of a hysteretic DC-SQUID has been 
measured at different values of the applied magnetic flux. At low temperature, the escape current 
and the width of the probability distribution are temperature independent but they depend on flux. 
Experimental results do not fit the usual one-dimensional (ID) Macroscopic Quantum Tunneling 
(MQT) law but are perfectly accounted for by the two-dimensional (2D) MQT behaviour as we 
propose here. Near zero flux, our data confirms the recent MQT observation in a DC-SQUID [3. 

PACS numbers: 



Quantum dynamics of a current-biased Josephson 
junction (JJ) have been extensively studied during the 
past two decades The phase difference across the J J 
obeys the same dynamics as the position of a fictitious 
particle moving in a one dimensional tilted washboard 
potential, whose average slope is proportional to the cur- 
rent bias. For bias below the critical current, the par- 
ticle is trapped in a local minimum. It can escape out 
by Thermal Activation (TA) or by Macroscopic Quan- 
tum Tunneling (MQT). Many experiments have demon- 
strated these two escape regimes. Quantized energy lev- 
els and resonant activation have been observed by ap- 
plying microwaves [J H| . Very recently Rabi oscillations 
and coherent temporal oscillations were observed in a 
current biased J J SJa] or in a Cooper pair transistor cou- 
pled to a biased Jjjfj. Experiments on RF-SQUID have 
also clearly demonstrated quantum dynamics in form of 
MQT pj, resonant tunneling or quantum superposition 
of states @ . 

By contrast and surprisingly, very few studies have 
been performed in a current biased DC-SQUID in the 
quantum regime. The dynamics of a DC-SQUID are de- 
scribed by two degrees of freedom and new phenomena 
can arise. The fictitious particle moves in this case in a 
two-dimensional (2D) potential. At zero magnetic field 
it has been predicted that escape out of a local minimum 
can occur via two different trajectories 0, @- At non 
zero magnetic field, MQT was derived only in the limit 
of small inductance of the loop 0, where the SQUID 
is nearly equivalent to a single junction with a flux- 
dependent critical current. In this case, the dynamics 
are one dimensional. To our knowledge only three previ- 
ous experiments were reported in the past discussing es- 
cape from a current-biased DC-SQUID. In the T A reg ime 
good agreement with 2D transition state theory [lfj has 
been obtained [Tl| . Similar results were obtained in an 
RF-SQUID I n the quantum regime no experimental 
evidence of two-dimensional behaviour has been demon- 
strated up to now. Sharifi et al. [13J claimed MQT ob- 
servations but neither the amplitude of quantum fluc- 
tuations nor the TA regime were understood. Recent 



measurements have clearly observed MQT and TA 
regimes in a DC-SQUID at zero magnetic field. 2D be- 
haviour could, however, not be observed because at zero 
flux the escape rate of their SQUID behaves exactly as 
that of a single Josephson junction, ie. a ID system. Un- 
derstanding the dynamics of a hysteretic DC-SQUID is 
very important in interpreting some recent quantum ex- 
periments flH ITfj| , or a recently proposed experiment of 
a one shot quantum measurement in a superconducting 
charge qubit using a DC- SQUID |l7j . 

In this letter, we report on an escape measurement in 
a DC-SQUID both in the thermal and in the quantum 
regime. In contrast to till now reported results on MQT, 
ID tunneling escape formula can not describe our exper- 
imental data. Since no theory has described MQT in a 
DC-SQUID in the 2D regime at non zero magnetic flux, 
we propose a method based on the zero magnetic field 
result derived in Ref. Q but taking into account the flux 
dependence of the 2D potential. Our experimental data 
support quantitatively this approach. 

The current-biased DC-SQUID consists of two Joseph- 
son junctions in parallel in an inductive superconducting 
loop. The whole system has two degrees of freedom de- 
scribed by ipi and ip 2 , the phase difference across the two 
junctions. Assuming identical JJs, each having a criti- 
cal current equal to Iq and a capacitance equal to Co, 
the dynamics of the system can be treated as that of a 
fictitious particle moving in a two dimensional potential 
given by 

U(x,y) = Uo (—sx — cosxcosy — rjsy + b(y — 7r<I>/$o) 2 ) 

(1) 

where Uq = ^qIq/it is the sum of the Josephson energy 
of the two JJs, s = I D c/(2Iq), x = (ipi + f2)/2, y = 
((pi — ip 2 )/2, b — $□/ (^ttLIq) and Iryc is the bias current. 
The asymmetry of inductances is parametrized by rj = 
(L-2 — L\)/L where L\ and L 2 are the inductances of the 
two branches and L = L± + L 2 is the inductance of the 
SQUID loop. $ is the applied flux and $o = h/2e is the 
flux quantum. 

The potential presents valleys and mountains with lo- 
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FIG. 1: Bidimensional potential of our SQUID which has pa- 
rameters s = 0.7, $/<I>o = -0.27 and b = 0.6 and 77 = -0.26. 
The particle and its escape path are drawn as illustration. 
Inset: photograph of the type of a DC-SQUID used in the 
experiment. 
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FIG. 2: Measured escape probability and escape rate (inset) 
as functions of the bias current at <E> / <E>o = 0.07, T — 42 mK 
and At — 50 /is. The dots are experimental data and the 
lines the MQT prediction of Eq. Values of I esc and A7 
are indicated by dashed lines. 



cal minima separated by saddle points along which the 
particle can escape. The barrier height between the min- 
ima and the saddle points depends on the bias current 
and magnetic flux and vanishes at the critical current I c . 
I c follows the usual flux dependence with <J> periodicity. 
In the following we consider b ~ 1 close to that of our 
DC-SQUID. In this limit, there exists only one optimal 
trajectory that passes through the saddle point. There- 
fore we will not discuss escape out of a local minimum 
through two different trajectories 0,0- Yet the 2D sig- 
nature still exists because the optimal trajectory is not a 
straight line along the x-direction but follows a path in 
the landscape with the tilt in the potential. 

In the limit s < s c where s c = I c /(2Iq) is the nor- 
malized effective critical current of the DC-SQUID, we 
can expand the potential (Eq. 0J) to the third or- 
der. The escape occurs in a current and flux depen- 
dent angle 9 with respect to the x-axis. Along this di- 
rection, we obtain the potential barrier height, AU = 
AUq\J s c cos 3 0/u(l + 77 tan 6 1 ) 3 / 2 , and the plasma fre- 
quency, uj p = w p o(MCOS0/s c ) 1 / 4 (l-|-r7tan0) 1 / 4 where u is 
the third derivative of the potential along the escape di- 
rection The quantities w p0 = 2 1 / 4 (7r/ c /$ C ) 1/2 (l- 
s/sj 1 / 4 and AU = (4^2/3) ($ /27r)J c (l - s/s c fl 2 are 
the usual plasma frequency and potential height, respec- 
tively, of a single JJ with an effective critical current I c 
and a capacitance 2 Co- The 2D character of the escape 
is exhibited in these expressions by the renormalization 
of the plasma frequency and the potential barrier height 
of a single JJ. Therefore the ratios lo p /lo p q and AU/AUo 
show deviations from ID behavior. 

Using the previous results of the renormalized plasma 
frequency and potential height, the escape rate out of a 
2D potential well, in the TA regime, T TA , is predicted to 
be Q 



where u> w ± and u> s ± are respectively the transverse fre- 
quency modes in the well and at the saddle point 
We numerically analyzed the modulation of the escape 
current as a function of flux in the TA regime taking 
into account the transverse frequency modes in the es- 
cape rate. For our SQUID parameters, the deviation of 
lu w j_/uj s j_ from unity has been estimated to be 2% at 
zero flux and it quickly decreases to zero when increas- 
ing I $/$o I- Since this small error appears only in the 
prefactor of Eq. , we will assume in the following that 
uj w ±/uj s ± = 1. 

At low temperature, the escape is dominated by quan- 
tum tunneling. This regime has been theoretically stud- 
ied at zero flux in 0, Q discussing the splitting of the 
two escape trajectories, or flux dependence in the low 
inductance limit (b >> 1) 9j. There is so far, to our 
knowledge, no theoretical treatment of flux dependence 
of 2D MQT, ie. in the regime b < 1. Therefore we pro- 
pose to use the MQT escape rate at zero magnetic flux 
pi Il4| in the limit b f» 1. We introduce in this escape 
rate, r M< 5 T , the 2D behaviour by taking the plasma fre- 
quency Up and the potential barrier height AU along the 
escape direction, analogously to what was done in the 
TA regime earlier [Tl| : 

t mqt = f 2D ^P./ 864tt AU/hujp exp (-36AU/5hv p ). 

(3) 

This 2D MQT rate thus depends on the magnetic flux 
and on the bias-current. In the following, because b « 
1, two escape trajectories are impossible and f2D — 1- 
Note, however, that the 2D behaviour remains through 
AU and u> p . At zero magnetic flux Y M Q T equals that of 
a single J J of critical current 2 Jo and capacitance 2Co. 

The escape probability P esc from the zero voltage state 
during time At is a function of the amplitude of the 
bias current, the flux and temperature. Measured escape 
probability using current DC-pulses of duration At = 50 
fis is plotted in Fig. |2| Each point corresponds to 4000 
current pulses. By increasing the bias current, P esc starts 
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from zero far below the critical current where the barrier 
is high, and it becomes unity at currents close to and 
above the critical one. In the following the escape cur- 
rent I esc is defined as the bias current corresponding to 
Pesc = 0.5 and the width of the escape probability dis- 
tribution AI as the difference between the currents at 
P esc = 0.9 and P esc = 0.1. The measured escape rate is 
extracted from T(Idc) = — ln(l — P eS c(lDc)) I 'At and it 
ranges from about 10 s _1 to 10 5 s _1 (inset of Fig. 0). 

The sample is a DC-SQUID consisting of two 
Al/AlO x /Al tunnel junctions obtained by e-beam lithog- 
raphy and shadow evaporation (inset of Fig. [TJ. The 
area of each J J is about 11 fim 2 . The tunnel resistance 
of the SQUID, Rt — 80 ft, and the superconducting gap 
A = 230 /xeV are deduced from the dissipative branches 
of the IV characteristic and yield Iq ~ 2.3 fiA using 
the Ambegaokar-Baratoff relation assuming two identi- 
cal JJs. From MQT of a single JJ with the same area, 
we deduced Co = 0.46 pF, which agrees with the typical 
specific capacitance of about 50 fF /(^m) 2 for an AlOx 
tunnel junction. The rectangular SQUID loop is about 
10 fim by 40 /im giving an estimated magnetic inductance 
of about 70 pH. The kinetic inductance of the SQUID, 
deduced from the thin film resistivity, p ~ 5.3/iOcm, is 
coarsely estimated to be about 120 pH. The total induc- 
tance L, the sum of the kinetic and magnetic inductances, 
is therefore about 190 pH. An inductance asymmetry 
exists because the two SQUID branches have different 
lengths. We estimate rj ~ —0.2 by assuming proportion- 
ality between length of a branch and inductance. The im- 
mediate electromagnetic environment of the DC-SQUID 
is determined by two 20 nm thin, 1 /xm wide and 700 
/im long superconducting aluminium wires terminated by 
two large pads for wire bonding. The estimated total in- 
ductance of the in-situ wires, L e = 4.6 nH, results again 
of the sum of 1.2 nH magnetic inductance of the copla- 
nar strips and of 3.4 nH kinetic inductance of the wires. 
Although these SQUID parameters have been indepen- 
dently estimated, the accuracy in determining them this 
way is not sufficient. Therefore we extracted the induc- 
tance of the DC-SQUID L = 244 pH and its asymmetry 
rj = —0.26, the zero-flux, the critical current Iq = 2.33 
fiA and the capacitance Co = 0.46 pF from the flux mod- 
ulation characteristics of the escape current in the low 
temperature regime (see Fig. |3J). These parameters are 
in very good agreement with the rough estimates given 
above. Our b = 0.6 is much smaller than b w 3 of the 
earlier works plU^ on MQT in a DC-SQUID, indicat- 
ing that the fictitious particle is much more affected by 
the two-dimentionality of the potential. 

The sample is enclosed in a cavity with resonances 
above 15 GHz, and it is anchored to the mixing chamber 
of the dilution refrigerator. This cavity is enclosed in a 
low temperature filter pp| consisting of thermocoax mi- 
crowave filters 0| and n-filters. The attenuation of the 
DC-lines is estimated to be at least 200 dB above 1 GHz 
at the mixing chamber temperature. Thermocoax filters 
were also introduced from 1.5 K to the low temperature 
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FIG. 3: Measured escape current and the width of the prob- 
ability distribution against external applied flux at low tem- 
perature (dots). ID MQT and 2D MQT prediction (Eq. 
are shown respectively by dashed and solid lines. 



filter, and low pass LC filters were used at room temper- 
ature. Current bias and voltage probe lines of the sample 
are separated from the computer using differential ampli- 
fiers. A superconducting shield was inserted inside the 
vacuum jacket of the cryostat and fj,- metal surrounds the 
dewar to protect the DC-SQUID from external flux noise. 
Using this measurement set-up, we separately measured 
a width of AI = 11 nA on a 1 /iA Josephson junction in 
its MQT regime. 

Figure [5] shows escape probability and escape rate at 
$/$ = 0.07 and T = 42 mK. At this particular flux 
value, the escape current is maximum (see Fig. [3J and the 
ratios u>p/u) p o and AU/AUq are equal to unity. Therefore 
the SQUID behaves as a single JJ. We can perfectly fit 
our data using 1D-MQT predictions (continuous line in 
Fig. J2J) without free parameters which demonstrates that 
1D-MQT theory describes the behaviour of a DC-SQUID 
at the maximum escape current. This first result con- 
firms the recent MQT observation at zero flux performed 
by Li et al, ^4|. It does not, however, demonstrate 2D 
escape behaviour. 

The escape current and the width versus applied ex- 
ternal flux are directly extracted from the measured 
switching current distribution and they are plotted in 
Fig. |3 in the MQT regime at T= 42mK . The width 
is periodic in flux with periodicity <i>o. It is minimum 
near <E» / <E»o = 0.07 and increases to reach maximum at 
<!>/<I>o — ±0.5. There is asymmetry between negative and 
positive flux values and discountinuities at &/&o = ±0.5. 
The usual ID tunneling model cannot account for the 
main feature of our results. Namely the ID model pre- 
dicts a decrease of the width, AI ~ (2/oSc) 3 ^ 5 m the 
MQT regime (dashed line in Fig. |3J which is contrary to 
our observation. Our proposed 2D formula in the MQT 
regime suggests increase of the width and it agrees with 
our data perfectly also quantitatively. It has no free pa- 
rameters since all of them were already fixed by the es- 
cape current characteristics. We claim that these results 
are the first experimental observation of 2D escape from 
the zero voltage state of a DC-SQUID in its MQT regime. 
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FIG. 4: Measured escape current and width of escape prob- 
ability distribution 81 vs. temperature (symbols) at $/$o = 
—0.25 together with predictions when including (continuous 
line) or excluding (dashed line) the influence of the environ- 
ment. We assume that the total escape rate is the sum of 
Tmqt and T T a- 

The escape current and width against temperature are 
shown in Fig. 0|for $/&o = —0.25. The escape current 
increases as temperature drops and saturates at low tem- 
perature. The width decreases with temperature and sat- 
urates at low temperature. These features are perfectly 
accounted for by Eq. © in the MQT regime and by 
Eq. J2J) in the TA regime (dashed lines). Yet the model 
seems to underestimate the measured width in the ther- 
mal regime by about 10 %, indicating that influence of 
the environment must be taken into account to better fit 
with the data. Finally, the observed cross-over tempera- 
ture (~ 80 mK) between MQT and TA is consistent with 
that predicted by the well-known formula htu p / 'finks) 

Q 



In order to estimate the effect of the environment on 
the escape process, we take into account the quadratic 
current fluctuations (SI 2 ). Since the SQUID is connected 
to the external circuit through inductive wires on the 
chip, we assume the electrical environment to be the 
inductance L e in series with an external resistor R ex t- 
For temperatures fc^T > HR ext / L e , the current fluctu- 
ations are given by (5I 2 ) = kBT/L e . For fc^T < haj p 
thermal fluctuations are slow as compared to the SQUID 
dynamics, and escape rate can be corrected by = 
T TA exp((<9A[/ /dI DC )(SI 2 ) /2{k B T) 2 ) in the TA regime 
and T A JQ T = T M( 3 T exp((d(AU/huj p )/dInn)(SI 2 )/2) in 
the MQT regime following the idea in [2l|. Using the 
estimated wires inductance, our data perfectly are fitted 
both in the quantum and in the TA regime by the fluc- 
tuation corrected expressions above (continuous line in 
Fig. 0} . We notice that like in the recent work in 01 the 
measured Al crosses over from MQT to TA much more 
gradually than what we predict. 

In summary, we have shown for the first time evidence 
of two dimensional macroscopic quantum escape of a hys- 
teretic DC-SQUID from its zero voltage state. Moreover, 
we have proposed a way to theoretically analyze the flux 
dependence of escape in the MQT regime, and we show 
that our data perfectly fit this model, also in the TA 
regime with the same values of the circuit parameters. 
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